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ABSTRACT 

Wc present a method to include the effects of gravity in the plasma physics code 
Cloudy. More precisely, a term is added to the desired gas pressure in order to 
enforce hydrostatic equilibrium, accounting for both the self-gravity of the gas and the 
presence of an optional external potential. As a test case, a plane-parallel model of the 
vertical structure of the Milky Way disk near the solar neighbourhood is considered. 
It is shown that the gravitational force determines the scale height of the disk, and 
it plays a critical role in setting its overall chemical composition. However, other 
variables, such as the shape of incident continuum and the intensity of the Galactic 
magnetic field, strongly affect the predicted structure. 

Key words: ISM:structure - methodsmumerical - gravitation - radiative transfer 
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1 INTRODUCTION 

Many astrophysical objects, from the scale of individual 
stars to the hot plasma of the intracluster medium, can be 
approximately described as being in hydrostatic equilibrium 
with the underlying gravitational potential. The density and 
temperature of the gas are related to the enclosed mass by 
the condition that the net acceleration vanishes. Thus, the 
hydrostatic equilibrium equation reads 



(1) 



1 dPjx) 
p{x) dx 

where P denotes the total (gas, magnetic, turbulent, cosmic- 
ray, and radiation) pressure, and x is the relevant spatial 
coordinate. In spherical coordinates, x = r, and the gravi- 
tational acceleration 



air) = 



GM(r) 



= ~— / x'^p{x) Ax 



(2) 



is determined by the total mass M(r) contained within ra- 
dius r. For a symmetric plane-parallel atmosphere, the rele- 
vant coordinate is the distance to the mid plane, z, and the 
gravitational acceleration is 



(3) 



g{z) = -27rGE(2) = -'i-wG / p{x) Ax 

Jo 

In either case, the solution to equation ([T| is 

P{x) = Po + APgrav(K) = Po+j p(u) g{u) Au 

Jo 



(4) 



where Pq is the central pressure. There are many situations 
of astrophysical interest where APgrav(a;) <^ Po, and the hy- 
drostatic equilibrium condition results in an almost isobaric 
atmosphere. However, depending on the total mass and the 
characteristic temperature of the gas, the gravitational field 
may be strong enough to have a significant effect on the 
structure of the system. 

One example of a photoionized gas where the gravi- 
tational acceleration is not negligible is the Galactic inter- 
stellar medium (ISM). Although interstellar matter only ac- 
counts for a relatively small fraction of the total disk mass, 
it plays a crucial role on the formation and evolution of the 
Galaxy. New stars are formed from the cold molecular phase 
of the ISM, and the energy input from the stellar popula- 
tion influences the structure and composition of the inter- 
stellar gas. The complex interplay between gas, stars, cosmic 
rays, and magnetic fields is responsible for the rich structure, 
dynamic s, and observa tional phenomenology of our Galaxy 
(see e.g. lFerriere|[200ll) . 

The present work describes the implementation of grav- 
itational hydrostatic equili brium in the plasma physics code 
Cloudy, last described bv lFerland et al.l (119981 ). Numerical 
details are given in Section [5] and Appendix |^ Section [3] 
discusses the results obtained for the vertical structure of 
the Milky Way disk in the solar neighbourhood, and a brief 
summary is given in Section 2] 



2 NUMERICAL IMPLEMENTATION 
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The Cloudy photoionization code computes the ionization 
structure, level populations, electron temperature, and ob- 
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servable spectrum of a gas, given its chemical composition 
and the spectrum of an ionizing source. The program allows 
the user to set the gas density at the illuminated face and 
to choose from several prescriptions to specify how it varies 
with distance. 

One of these prescriptions is constant pressure, which 
may refer to the gas or the total pressure 

P{x^ = -Pgas ~t~ -Pram "4" -fturb + -Pmag + -Plincs + ^-Prad (5) 

where Pgas = nkT is the thermal pressure, -Tram — P^wind 
and Pturb = arise from uniform and turbulent 

motion, respectively, Pmag = is the magnetic pres- 
sure, Piines is the radiation pressure of trapped emi ssion 
lines l|Ferland fc ElitzuJlTgsi : lElitzur fc Ferlandlligsd ) and 
dx compensates for the acceleration Orad 
due to the absorption of the incident radiation. This term is 
analogous to APgrav and, to some extent, one could argue 
that the command should already be named "hydrostatic" 
rather than "constant pressure" even if it did not account for 
the effect of gravity. Since Orad is positive (the momentum 
of the absorbed photons points outwards), AP^ad induces 
a positive pressure gradient (see Figure [3|, whereas gravity 
must be balanced by a negative gradient. 

We have modified the constant pressure option in 
version 08.00 of Cloudy to account for the gravitational 
term APgrav. Source code is available upon request, and a 
brief user guide is provided in Appendix [X] The required 
modifications are relatively minor, and they are not expected 
to incur in any side effect on the functionality of the algo- 
rithm. They are currently being implemented in the devel- 
opment version of the code and will be available in the next 
official releas^. 

First, new variables are created to store the values of 
APgrav (initialized to at the illuminated face), the sym- 
metry of the problem (spherical or plane-parallel) and the 
mass (0 by default) of an additional component located at 
the centre. 

As Cloudy moves away from the illuminated face. 



APgrav is updated according to 
APi+i = APi + pi Qi hi 



(6) 



where APi and APi+i denote the old and new values, while 
Pi, Qi, and hi are the density, gravitational acceleration, and 
thickness of the last zone, respectively. The gravitational 
acceleration is computed according to expression ((2]) or ((3]) 
depending on the geometry. Cloudy keeps track of the cu- 
mulative mass M{r) and the gas column density S(z) in each 
case. Additional mass components or any external (gravita- 
tional or not) potential can be trivially accounted for, simply 
adding the appropriate contribution to the local acceleration 

Finally, one has to modify the pressure loop so that the 
the term APgrav is added right after APrad to balance the 
gravitational force. Then, the code proceeds exactly as in 
the original constant pressure case, trying different values of 
the gas density until the deviation from the desired pressure 
is smaller than the requested tolerance. 



1000 \r 
100 
10 \r 
1 

0-1 \r 
0-01 

10"- 

10"' 
1000 \r 

100 

10 \r 

1 

0.1 \r 

0.01 
' 10"^ 
10"' 



" 1000 
C 100 t- 



0.1 

0.01 ir 
10"- 
10"' 
1000 
100 



0.1 

0-01 \r 

10 
10"' f 



n-, = 1.0 cm 



Hn = 0-5 cm 



□ n = 0.1 cm 



10 100 1000 

z [pc] 

Figure 1. Molecular, neutral, and ionized Hydrogen density as 
a function of heiglit z with respect to the disk mid plane (dotted, 
solid, and dashed lines, respectively) for different values of the 
mid- plane density no. 



3 THE MILKY WAY DISK 

In order to test the proposed implementation of hydrostatic 
equilibrium, we consider a simple model of the ISM in the 
solar neighbourhood, assuming a plane-parallel geometry 
that is symmetric with respect to the G alactic mid plane, 
z = 0. The gas chemical composition (ICowie fc Songailal 
ll986l : ISavage fc SembachI 19961: iMever et al.lll998l ). including 
dust grains ( Mathis et al.lll977 ). is set by the abundances 
ISM command, w hile cosmic ray heating and i onization are 
treated following iFerland fc Mushotzkvl ([l98J) with the in- 
struction cosmic ray background (see HazjQ, the Cloudy 
documentation, for details). 

Given a mid-plane density no, the vertical structure 
of the disk is self-consistently computed by the program 
under the assumptions of ionization, thermal, and hydro- 
static equilibrium. An important limitation is that the gas 
is described by a single density and temperature that may 
vary with distance, but the coexistence of different phases at 
the same point is not implemented (see e.g. |Boul a. res fc Cm3 
ll990l : ISilich~fc" Tenorio- Taglelll998l : iFerrierelliooTI : [cQi2005l . 



^ The current stable version of Cloudy can be downloaded from 



|http : / /wiki . nublado. org / wiki /Download | 



Available at |http:/ /www. nublado. or"g| 
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Figure 2. Gas temperature (top) and pressure (bottom) along 
the 2;-axis for values of the mid-plane density no = 1.0, 0.5, 0.2, 
and 0.1 cm~^. 



for examples of multiphase models of galaxy disks in hydro- 
static equilibrium). Another simplification is that the radi- 
ation field is assumed to strike the system from one side 
(in our case, the mid plane), whereas in reality the stars 
are distributed across a scale height comparab le to that of 
the g aseous component (e.g. Miller fc Cox. 19931 ; iRobin et al.l 
I2OO3I ). With these caveats in mind, it is not surprising that 
our simple model does not provide a perfect match to the 
available observational data, ft should be able, though, to 
give a rough estimate of the volume-averaged densities and 
scale heights of the different ISM phases. More importantly, 
this kind of models offers some quantitative insight on the 
role played by the different physical mechanisms that may 
be at stake. Although the present study is mostly concerned 
with gravity, we will also explore the influence of the spec- 
tral energy distribution of the incident radiation, as well as 
the presence of an external potential and a magnetic field. 



3.1 The role of gravity 

To illustrate the effect of gravitational acceleration on the 
equilibrium structure of a generic photoionized region, we 
ran a first model where the incident continuum consists of 
the cosmic microwave background and the interstell ar ra- 
diatio n field, interpolated from the measurements of iBlackl 
1 19871 ) by means of the table ISM command. 

Results for different values of the central density no are 
plotted in Figure [T] In all cases, there is an inner region 
where the gas is almost fully ionized, a thin layer of molec- 
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Figure 3. Density of the different Hydrogen phases (top), elec- 
tron temperature (middle) and total pressure (bottom) in the 
standard constant pressure case. 



ular and atomic Hydrogen at high densities, and an outer 
region of low-density gas where Hydrogen is predominantly 
neutral, becoming more and more ionized at large distances 
as density decreases. Gas temperature and pressure as a 
function of z are shown in Figure [51 The temperature of 
the cold dense phase is of the order of 10 — 100 K, while 
the diffuse gas ranges from 5000 to 10000 K. Gas pressure 
is roughly constant up to the ionization front, but it de- 
creases rapidly through the molecular layer due to the high 
gas densities (above 100 cm~^). The pressure stabilizes again 
in the outer region, although the gravitational acceleration 
imposes a negative pressure (and density) gradient. The ef- 
fect increases with distance; at some point, the cumulative 
gas mass and therefore the gravitational force are so high 
that both density and pressure decay extremely fast. Gas 
temperature, on the other hand, increases steadily towards 
the Galactic halo. 

This structure does not provide a realistic description 
of the local ISM. On the one hand, the approximation that 
all the radiation comes from an infinitely thin disk at the 
mid plane creates a photoionized region at low z. Although 
there is indeed a thin layer of ionized gas near the Galactic 
plane, c omposed of both disc rete Hll regions and diffuse gas 
(see e.g. iPaladini et al]|2005l ). its contribution near the so- 
lar neighbourhood is negligible (e.g. iTavlor fc Cordeslll99^ : 
ICordes fc Lazidl2002l ). On the other hand, neither the po- 
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Figure 4. Incident continuum given by the table ISM command 
(dashed Une), absorption of ionizing photons by the extinguish 
command with a Hydrogen column density of 10'^^ cm~^ (solid 
line), and propagation through a uniform gas slab with constant 
density of 1 cm~'^ and the same Hydrogen column density (dotted 
line). 

sition nor the thickness of the neut ral phases afiree with 
observations of the Milky Way ISM jBoulares fc Coxl Il990l : 
lFerriere|[200ll : iNakanishi fc Sofuell2003l . l2006l ). Nevertheless, 
a central ionized region may be present in several astro- 
physical scenario^, and therefore it is interesting to test the 
influence of gravity in this particular example. 

The role of gravity can be assessed by comparison with 
the standard constant pressure solution, depicted in Fig- 
ure [3] for a mid-plane density no = 0.5 cm~'^. The vertical 
structure of the disk is basically the same (i.e. gravity does 
not play a significant role) up to the ionization front. There, 
the temperature decreases sharply, and the density has to 
rise accordingly to compensate. In the standard implementa- 
tion, pressure increases in order to balance the radiative ac- 
celeration (APrad > 0). The effect is very significant near the 
ionization front, but it rapidly decreases as photons become 
absorbed. Although the pressure keeps increasing monoton- 
ically with z, the slope is quite shallow, and the physical 
state of the gas in the outer region (density, chemical com- 
position, and temperature) remains approximately constant. 
Thus, the main effect of gravity on a plane-parallel atmo- 
sphere is to set the extent of the molecular layer and the to- 
tal scale height of the system. Without gravity, the neutral 
(mostly molecular) layer would simply extend to infinity, as 
can be readily seen in Figure O 

3.2 The incident radiation field 

As discussed in the documentation of the table ISM com- 
mand, ionizing radiation between 1 and 4 Rydberg is heavily 

^ Note that such an equilibrium configuration is, in any case, 
Rayleigh- Taylor unstable. Cold gas blobs and filamentary struc- 
tures would slowly detach from the dense, neutral material and 
fall towards the low-density, ionized region, whereas hot gas bub- 
bles would tend to rise buoyantly through the neutral phase. 
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Figure 5. Effects of the incident continuum on the vertical struc- 
ture of a model with no = 1 cm""^. The table ISM case is depicted 
in Figures[l]and[2l The two top panels show the molecular, neutral 
and ionized Hydrogen densities for the incident radiation fields 
obtained by applying the extinguish command and the trans- 
mission through a uniform slab. The third and fourth panels show 
the gas temperature and the total pressure in each case. 

absorbed by the neutral Hydrogen in the ISM so that few 
such photons exist, at least in the Galactic plane. In order 
to describe the actual radiation field incident on a typical 
region in the Galactic plane, the use of the extinguish com- 
mand is recommended. However, it is also stated that this 
command should not be used except as a quick test, and a 
more physical extinction of the continuum may be accom- 
plished by transmitting it through a model of the absorbing 
slab. 

For the sake of simplicity, we assume a uniform density 
of 1 cm~'^ and a total column densit y of 10^^ cm~^, appro - 
priate for the Milky Way disk (jPickev fc LockmarJ il990l ). 
The main difference between the resulting continua, shown 
in Figure |4l is that the extinguish command does not affect 
the energy band between 0.1 and 1 Rydberg, responsible for 
molecular Hydrogen photodissociation. When the incident 
radiation field is propagated through the ISM, many of these 
ultraviolet photons become absorbed by the dust component 
and re-emitted in the far infrared. 

The details of the interstellar radiation field have a dra- 
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Figure 6. Total stellar mass density (solid line) and contribu- 
tions of indiv idual populat i ons w ith diflferent ages (dotted lines) 
according to iRobin et alj l l2003tl . The approximation given in 
equation (0 is shown as a dashed line. 



matic impact on the vertical structure computed by the 
code. Figure [5] compares the results obtained for no = 
1 cm~'^ using the extinguish command and those obtained 
by propagating the table ISM field through the uniform slab 
before reaching the Galactic mid plane. The main difference 
between either method and the unextinguished case repre- 
sented on the top panel of Figure[T]is the absence of the inner 
layer of ionized gas. Central temperatures and pressures are 
almost two orders of magnitude lower, and, since there is no 
phase transition associated to any ionization front, the den- 
sity of the cold neutral material is equal to no, also orders 
of magnitude lower than in the model without attenuation. 

The absorption of the ultraviolet photodissociating ra- 
diation turns out to be very important as well. The location 
of the transition from the cold to the warm neutral phase 
may change by a large factor, but the most relevant issue 
is the predicted fraction of Ifydrogen in molecular form. In 
particular, the extinguish command overestimates the pho- 
todissociation rate so much that there is virtually no molec- 
ular Hydrogen at all. Our results can only strengthen the 
advice given in Hazy of modelling the absorbing material 
rather than using the extinguish command. 



3.3 Stellar mass and Galactic magnetic field 

We have shown that photoionized regions must have a fi- 
nite extent because of self-gravity. The exact value of the 
characteristic scale length depends, however, on many other 
factors. For instance, the gas may represent only a small 
fraction of the total mass of the system, and there can 
be a significant amount of non-thermal pressure support 
from turbulent motions, cosmic rays, and magne tic fields 
l|Abel fc Ferlandll2006l : iPeUegrini et ai]|2007l . |2009l) . 

In the case of the Milky Way, the distribution of stel- 
lar mass in the Galactic disk can be inferred by combining 
number counts of individual stars and dy namical arguments . 
We approximate the results obtained bv lRobin et al.l (|2003l ) 
with the expression 
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Figure 7. Densities of the different Hydrogen phases (top pan- 
els), temperatures (third panel) and pressures (bottom panel) for 
two models including a stellar mass component of the form (0. 
In addition, one of them also includes a tangled magnetic field 
with an intensity Bg = 5 /iG at the Galactic plane. 



p. (2) = Po 



0.9 



2z 



z 
hi 



(7) 



where po = 0.042 M0 pc~'^ and /i* = 520 pc. This function, 
plotted as a dashed line in Figure [S] provides a reasonable 
approximation for z < h,, well beyond the half- width half- 
maximum of the distributiorQ. At larger distances, we take 
p, ~ 0. The adopted value of the central stellar density 
Po implies the same mass as 1.7 cm"'' of pure Hydrogen, or 
~ 1.2 cm~^ in the case of the ISM, where Hydrogen accounts 
for ~ 72 per cent of the gas mass. The contribution of stars 
to the local density near the Sun is thus comparable to that 
of the gaseous component. The cumulative surface density 



E*(z) = 2po 



(8) 



increases up to ~ 17.5 M© pc at 2: = /i, , about 3 times 
the observed surface density of Hydrogen (in all forms) at 



^ In fact, p«(h«) = O.lpo 
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the solar radius, and 2 times that of the ISM gas, adding up 
the mass in Hehum and heavier elements. 

When the contribution of the stellar component to the 
gravitational acceleration is taken into account, the pres- 
sure gradient required to maintain hydrostatic equilibrium 
increases, and the transition between the cold and the hot 
phase moves inwards by about a factor of 2 (top panel in 
Figure[7ll. In general terms, the presence of an external grav- 
itational potential will always push the gas distribution in- 
wards, leading to a reduction of the relevant scale lengths 
and the total gas mass obtained for the system. 

On the other hand, magnetic fields tend to counter- 
act the action of gravity. Due to the effect of line dragging, 
the intensity of a tangled field increases with gas density as 
B P C v?^^ (observationally, B oc n"', w ith 0.5 < k < 1, see 
e.g. ICrutchej|T999l : iHennev et al]l2005l ). and the associated 
magnetic pressure varies as Pmag oc B^ oc n*^^. This equa- 
tion of state is stiffer than that of an ideal gas, and therefore 
a much milder change in gas density is necessary in order to 
achieve hydrostatic equilibrium. 

The second panel of Figure[7]shows how the inclusion of 
a tangled magnetic field with mid-pla ne intensity Bo — 5 
l|Beck[|200ll : iHeiles fc Crutchedl20ol ) changes the structure 
of our ISM model. As can be inferred from the bottom pan- 
els, the magnetic field provides most of the pressure sup- 
port against gravity. The extent of the cold phase increases 
by roughly an order of magnitude, reaching now hundreds 
of pc. In the outer, hot layer, pressure is dominated by the 
thermal component, and the magnetic field (for our adopted 
value of Bo) does not play a significant role. 



3.4 Comparison with observations 

Our single-phase model can only describe the average prop- 
erties of the ISM. More precisely, it should be able to predict 
the density of molecular, atomic, and ionized Hydrogen as a 
function of the distance z from the mid plane of the Galaxy, 
averaged over some scale larger than the typical size of the 
inhomogeneities . 

Much observational work has been devoted to the deter- 
mination of the average vertical structure of the Milky Way 
ISM near the solar neighbourhood. For the ionize d compo- 
nent, we will adopt the distribution derived by iRevnolclsl 
lll99iri 



= 0.015 exp 

cm~^ \ 



70 pc 



1+0.0025 exp( 



900 pc 



(9) 



from pulsar dispersion measures, where the first term rep- 
resents the contribution from discrete Hll regions near the 
Galactic plane, and the second term corresponds to the dif- 
fuse ionized gas. For the neutral Hydrogen, we use the for- 
mula 



nHi,H2(-2=) = '^0 sech 



Zl/2 



■ ln(l + V2) 



(10) 



with no — 0.6 cm~'^ and 2:1/2 = 12 5 pc for the 
atomic component (|Nakanishi fc Sofuell2003l ), whereas no — 
0.09 cm~^ and Zi/2 = 93 pc describe the molecular gas 
l|Nakanishi fc Sofud 20061 ). The corresponding density distri- 
butions are shown in Figure|8]by dashed lines, and the effect 
of 20 per cent errors in each normalization and scale length 




z [pc] 

Figure 8. Comparison between the model predictions and the 
observed densities of molecular (top), atomic (middle) and ion- 
ized (bottom) Hydrogen. The observational data are represented 
as dashed lines, with the grey areas delimited by dotted lines cor- 
responding to 20 per cent error bands. Black solid lines display 
the results of the model with no = 1 cm~'^, Bo = 5 fiG, and 
stellar density jTj- A model with no = 0.75 cm~^. Bo = 4 ^G, 
K = 0.5, and N(R) = 10^^ cm~^ is plotted in solid blue (see text 
for details). 



is represented by shaded areas enclosed by dotted lines. This 
estimate roughl y corresponds to the actual errors quoted in 



Reynolds (1991), and it is substantially smaller than those in 
Nakanishi fc Sofue (.2006. ). No information about the obser- 



vational errors is provided in iNakanishi fc Sofud l|2003h . but, 
given the large dispersion in t he Hydrogen column dens ities 
along different lines of sight (|Dickev fc Lockmarj|l990l ). as 
well as the amplitude of the fiuctuations in the azimuthally- 
averaged surface density profiles, 20 per cent is arguably a 
realistic estimate. 

The results of a model with central density no — 
1 cm~^, central magnetic field Bo = 5 /iG, and a stellar 
distribution given by equation (O are also shown in Fig- 
ure [8] as black solid lines. Although the model gives a rea- 
sonable order-of-magnitude estimate (in agreement with the 
observed values within a factor of two) of the scale height of 
the Galactic disk, as well as the densities and total masses of 
both atomic and molecular Hydrogen, there are some impor- 
tant discrepancies: first, the molecular fraction near the mid 
plane is severely underestimated; second, the central density 
is almost a factor of 2 above the observed one; and third, 
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the amount of diffuse ionized gas is extremely low over the 
whole radial range, especially in the inner regions. 

The first two problems are closely related. One could 
easily specify a lower mid plane density, but then the amount 
of molecular Hydrogen would decrease to unacceptably low 
levels. It seems that, even after transmission through a col- 
umn density of 10^^ cm~^, the intensity of the photodis- 
sociating continuum is still too high for molecules to form 
at the observed rate. The third problem is also related to 
the simple prescription adopted for the incident continuum, 
since almost all the ionizing radiation becomes absorbed be- 
fore reaching the illuminated face. In reality, young stars are 
spread around the Galactic plane. Adding any thin source 
of ionizing photons at the mid plane results in an Hll re- 
gion qualitatively similar to the one observed in Figure [T] 
in order to reproduce the observed ionization fraction, one 
would need to specify a distributed radiation source, but 
such feature is unfortunately not currently implemented in 
Cloudy. 

The values of the model parameters can be tuned to 
better match the observed total density profiles and the 
molecular Hydrogen fraction. Blue solid lines in Figure [S] 
shows the results of a model where the mid-plane density is 
set to no = 0.75 cm~^, the intensity of the magnetic field 
is _Bo = 4 at the centre and varies with gas density as 
(B/Bo) ~ (n/no)°'^, and the incident continuum has been 
propagated through a column density of 7V(H) — 10^^ cm~^. 
Nevertheless, we strongly advise caution when interpreting 
these results. Although the fair agreement obtained for the 
atomic and molecular data provides encouraging support for 
the validity of our implementation of gravity, a more real- 
istic model of the ISM is clearly required in order to derive 
robust quantitative constraints on its physical properties. 
Based on our crude model, we would simply argue that rea- 
sonable values of the input parameters (in particular, central 
gas density and magnetic field) yield reasonable scale heights 
for the disk. Atomic and molecular Hydrogen densities typ- 
ically agree with observations within a factor of two, but a 
prescription for distributed sources is necessary in order to 
model the ionized component. 



4 CONCLUSIONS 

The present work describes an implementation of gravity 
in the well-tested photoionization code Cloudy. For an at- 
mosphere in hydrostatic equilibrium, the gravitational force 
has to be balanced by a pressure gradient according to equa- 
tion (in. 

By way of example, we consider a plane-parallel model 
of the Galactic ISM in the solar neighbourhood to illustrate 
the infiuence of the gravitational acceleration on the physical 
properties of the gas. The main result, valid in the general 
case, is that gravity sets the scale height of the system. If 
self- gravity is not accounted for, there is nothing that pre- 
vents the gas (mostly molecular at large radii) from extend- 
ing to infinity. When gravity is present, the pressure has to 
decrease in order to maintain hydrostatic equilibrium, and 
the gas becomes less dense. The precise structure depends 
on many other variables; we have investigated the effect of 
the incident radiation field, the presence of an external mass, 
and a tangled magnetic field, but other ingredients, such as 



turbulent motions, cosmic rays, dust grains, or variat ions 
in the atomic rate coefficients (see e.g. Abel et aLll2008l ), as 
well as the external pressure of the intergalactic medium 
( Sil ich fc Tenorio-Tagld ^0^) may also play an important 
role. To summarize, the present results show that gravity 
(in particular, the self-gravity of the gas) should not be ne- 
glected in studies of the structure of the atomic and molec- 
ular layers of a photoionized region, especially when one is 
interested in their total extent. 

Concerning the specific problem of the vertical structure 
of the Milky Way ISM, the models discussed here are able to 
reproduce the observed distribution of atomic and molecular 
Hydrogen with an accuracy of about a factor of two, but the 
amount of diffuse ionized gas is severely underestimated at 
all heights. Much more realistic results could be obtained by 
providing support for distributed radiation sources, as well 
as a treatment for the coexistence of different gas phases at 
the same point. 
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APPENDIX A: USER GUIDE 

In order to specify a hydrostatic equilibrium configuration, 
the constant pressure command must be included in the 
Cloudy input script. In addition, a gravity command with 
syntax 

gravity <symmetry> [<factor> [_log]] 

has been implemented. The compulsory argument symmetry 
must be either spherical or plane-parallel in order to tell the 
code whether to use equation ((2} or (|3} to compute the grav- 
itational acceleratiorjj. An optional factor can be introduced 
that multiplies the gas mass, mimicking the presence of an 
external component with identical distribution. A more flex- 
ible way to specify an external potential is provided by the 
command 

gravity external <mass> [<extent> <index>] [_log] 

that models a distribution of mass around the centre of sym- 
metry of the system. In the spherical case, the first number 
represents the mass (in solar masses) of a pointlike object 
located at r = 0; in the plane-parallel case, it corresponds 
to the total surface density (in Mq pc~^) of a uniform thin 
sheet at the mid plane 2 = 0. In both cases, the total mass 
mo is assumed to be distributed over an extent xq as a 
power law with index a, i.e. m(x) — mo{x/xo)". All num- 
bers are interpreted as linear, unless the keyword -log is 
specified. 

^ Actually, the command parser looks for tiie character sequences 
"sphe" and "plan", so expressions like "sphere" or "planar" are 
also valid. 



The following is an example of a typical input script 
used in the preparation of this work: 

# 

# Cloudy input file for hydrostatic equilibrium 

# Yago Ascasibar (UAM, Fall 2009) 
# 

# 

c 

title - Hydrostatic equilibrium - 

c 

c Incident continuum 

c 

# CMB 

# table ism 

# extinguish column = 21 leak = 
table read "attenuated_ism.dat" 
nuF(nu) = -2.0453234 at 7.541e-05 Ryd 

c 

c Density and chemical composition 

c 

hden 1 linear 
abundances ism 

c 

c Gravity 

c 

constant pressure 
gravity plane-parallel 

gravity external Sigma = 43.680, zO = 520, a = 1 
gravity external Sigma =-39.312, zO = 520, a = 2 
gravity external Sigma = 13.104, zO = 520, a = 3 

c 

c Non-thermal components 

c 

cosmic ray background 
magnetic field -5.3 

c 

c Stopping conditions 

c 

stop temperature off 

stop thickness le4 linear parsec 

stop eden -4.5 

iterate to convergence 

c 

c Output files 

c 

punch overview "overview.txt" last 
punch pressure "pressure.txt" last 

# punch transmitted continuum "trems.txt" last 

# punch hydrogen conditions "Hydrogen.txt" last 

c 

c ... Paranoy® Ruiz! ; ~D 
c 

The first blocks set a label to identify the run, the in- 
cident radiation field, the central Hydrogen density no and 
the gas chemical composition. Commands related to hydro- 
static equilibrium come next. Note, in particular, how sev- 
eral gravity external commands have been combined in 
order to specify the stellar mass distribution given by expres- 
sion The integration continues until one of the stopping 
conditions is met (the disk thickness reaching 10 kpc or the 
electron density falling below ~ 3 x 10~^ cm~^), and it is 



iterated until the optical depths converge. Finally, the main 
physical properties of the gas (e.g. densities, temperatures) 
as a function of distance, as well as the different contribu- 
tions to the total pressure, are printed to output files on the 
disk. 



Self- gravitating photoionized gas 9 



